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We show that a constriction-type Josephson junction realized by an epitactic thin film of a d-wave 
superconductor with an appropriate boundary geometry exhibits intrinsic phase differences between 
and 7r depending on geometric parameters and temperature. Based on microscopic Eilenberger 
theory, we provide a general derivation of the relation between the change of the free energy of the 
junction and the current-phase relation. From the change of the free energy, we calculate phase 
diagrams and discuss transitions driven by geometric parameters and temperature. 
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I. INTRODUCTION 

The change of the free energy of a Josephson junction 
(JJ) evoked by the variation of the phase determines the 
intrinsic phase difference in the unbiased ground state. 
Usually, the coupling energy between the electrodes of a 
JJ is positive and the current-phase relation is sinusoidal, 
corresponding to a vanishing intrinsic phase difference. 
For the peculiar case of negative coupling, however, in- 
trinsic phase differences of ir are possible (see Refs. [lIQ 
and references therein). 

In the crossover regime between positive and nega- 
tive coupling, higher harmonics dominate the current- 
phase relation. This behavior has been studied for 
different types of Josephson devices: tunneling and 
in particular grain boundary JJs involving d-wave su- 
perconductors^^, controllable superconductor-normal 
metal-superconductor JJs (SNS) 6 , superconductor- 
ferromagnet-superconductor JJs££, periodically alter- 
nating 0-7T JJs2*i£ and grain boundary JJs in noncen- 
trosymmetric superconductors^. The dominating higher 
harmonics close to the 0-7T crossover lead to additional ze- 
ros of the current-phase relation. Whether or not these 
additional zeros are related to stable energetical minima, 
and accordingly to intrinsic phase differences of neither 
nor 7r, can only be decided by consideration of the free 
energy of the J J. 

In the present work, we show that intrinsic phase dif- 
ferences in the full range < 70 < tt occur across d- 
wave superconductor microbridges. The geometry un- 
der consideration consists of a stripe of a c-axis oriented 
epitactic thin film of a d-wave superconductor, which is 
narrowed down from one side by a wedge-shaped inci- 
sion (Fig. [1]). Accordingly, the Josephson effect follows 
solely from the lateral constriction of the thin film and 
emerges if the width w of the bridge is of the order of 
the coherence length £0 of the superconducting material 
or below. Such a microbridge configuration is clearly 
distinct from a grain boundary, superconductor-isolator- 
superconductor or SNS tunneling J J 1213 . It should be 
emphasized that in the geometry under consideration, 
Fig. [H there exists no grain boundary. According to the 
terminology introduced in Ref. [14], this microbridge con- 




FIG. 1: (Color online) Geometry denning the microbridge 
Josephson junction based on a c-axis oriented epitactic thin 
film of a d-wave superconductor. A stripe of the supercon- 
ducting material is narrowed down from one side by a wedge- 
shaped incision with the opening angle f3. The width of the 
resulting microbridge-like junction is given by w. Two typ- 
ical quasiparticle trajectories across the junction: 1-without 
reflection, 2-with reflection. 

figuration belongs to the weak link type JJs since the 
electrodes are not electrically separated by a tunneling 
barrier. 

In order to characterize the dc Josephson effect, we 
calculate current-phase relations based on microscopic 
Eilenberger theory. For the derivation of the intrinsic 
phase difference, we make use of the relation between 
the change of the free energy of the junction £(7) — £(0) 
and the current-phase relation 1(7). We give a very gen- 
eral derivation of this relation which is valid for arbitrary 
structures exhibiting a current-phase relation in the full 
temperature range < T < T c as well as in the pres- 
ence of an external magnetic field. The intrinsic phase 
difference of the microbridge will be discussed in terms 
of phase diagrams, justified by the thoroughly derived 
change of the free energy. 

II. BASIC MECHANISM 

The intrinsic phase shift of the device is a direct con- 
sequence of the d-wave symmetry. If the width of the 
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junction w is large, quasiparticle trajectories without and 
with a reflection at the straight edge opposite to the 
wedge contribute to the total current across the junc- 
tion (trajectories of type 1 and 2 in Fig. [1]). If the 
constriction is narrow enough, however, the dominant 
contribution to the total current stems from trajecto- 
ries which get reflected (type 2). If the orientation of 
the d-wave is for example a = 7r/4 (nodal surface), all 
reflected trajectories suffer a sign change of the pairing 
potential which leads to the formation of pronounced zero 
energy Andreev bound states at the surfac o 15 i 16 i 17 . For 
< a < 7r/4, still a fraction of all trajectories contributes 
to the formation of zero energy Andreev bound states, 
which engender anomalous counterflowing quasiparticle 
surface current o 18 i 19 and also intrinsic phase shifts^. 

III. THEORY 

A. Current-Phase Relation 

In order to calculate current-phase relations of the J J, 
we employ microscopic Eilenberger theory 2 ^. The self- 
consistency equation which has to be solved for the pair- 
ing potential A(r, kp) reads 

/ d2k F l V P™r} k F ,k' ^ 2(2 /i\ 

A(r,kp)= / — * — —. - 2nk B T > (1) 

V J J FS (27T) 3 \hv' F \ eZ^0 1 + ab 



Here, FS is the Fermi surface, = Vi?(k' F ) is the Fermi 
velocity, [V pa i r ]k F ,k' F is the pairing interaction matrix, 
e n = {2n J rl)irkBT are fermionic Matsubara frequencies, 
and a = a(r, k^,£ n /) and b = &(r, k. f F ,e n ') are the Ric- 
cati amplitude o 22 ! 23 . The selfconsistency equation for the 
current density j(r) is given by 

u r) - 2e f ^ tl^T ^ , l^ab 

i{r) - 2e L(^r 1^1 3^ 

The selfconsistency equations ([1]), ([2]) allow for the mi- 
croscopic calculation of the current-phase relatio n 20 ! 24 . 
Selfconsistent solutions guarantee current conservation 
(V • j(r) =0), but can in general only be found numeri- 
cally. 



B. Free Energy 

In order to derive the relation between the change 
of the free energy of the junction £ (7) — £ (0) and the 
current-phase relation 1(7), we start from the Eilen- 
berger functional 21 for the free energy £ (A, A^, A; a, 6), 
parametrized by the pairing potentials A, A^, the vector 
potential A(r) and the Riccati amplitudes a, b: 
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Here, A = A(r') is the total magnetic vector poten- 
tial, T$ext = B e xt(r / ) is the external magnetic field, and 
a = a(r', k' F , e n ') and b = 6(r', k' F , e n >) are the Riccati 
amplitudes. 

Consider a general variation of this functional: 



BE 

d£ = 8A dA <8* 
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The variation with respect to A^ yields the selfconsis- 
tency equation for A(r, kp), Eq. (pQ), and the variation 
with respect to A a corresponding selfconsistency equa- 
tion for At(r, kp): 



V F) J FS (2tt) 3 \hw' F \ B ^ n l + ab K ' 



The variation with respect to b yields the Riccati differ- 
ential equation for a and vice vers a 22 i 23 : 
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Accordingly, all variations vanish in the case of a self- 
consistent solution, i.e. at the stationary point of the 
functional, independent of the gauge. 
After making use of the selfconsistency equation for the 
currents, Eq. ([2j), and identification of the external cur- 
rents via rotBe^t = fiojext, we find for d£/dA: 



d£ = I d 3 r' 



(j(r') 
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If one integrates over all space, the variation of £ 
with respect to A vanishes due to Maxwell's equation, 
rot rot A = /i (j + iext)- Restricting the integration to 
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a finite volume V enclosing the junction, the external 
currents ] ext drop out. Now we carry out the gauge 
transformation {A -> Ae^, A 1 " -> A+e - ^, a -> ae^, b -> 
&e"^, A->A+ Ay0}, leading to 



- / 



dV 



j(r') - — rot rot A(r / ) 1 d(V0) (7) 
Mo J 



with the gauge-invariant phase <\> = (j) + ^ dl • A. 
Considering the volume V according to Fig. [2] and using 
basic vector calculus, it can be shown that the contri- 
bution of the vector potential vanishes. Integration by 
parts, exploitation of current conservation and applica- 
tion of Gauss's theorem then results in 



d£ = 



h_ 

Ye 



da n • j(r') dcf) 



(8) 



S=dV 



Only the parts of the surface S = dV where the current 
enters into or leaves the volume V contribute. Since cross 
sections of the superconductor far from the junction are 
planes of constant gauge-invariant phase, 
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with the total current / and the gauge-invariant phase 
difference ^ = 4>r — <\>l — j[ J** d\ • A and finally 



£(7) 



m = Ye£ dl ' nj,) 
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Because of current conservation, the total current for the 
current-phase relation can be taken at any cross 

section of the superconductor S sc : 



7( 7 ) = j ^'n-j(r') 

Js sc Jfs(^) 3 IK- 1 ^ 



(11) 

1 — ab 



Quasiparticle bound states as well as the supercurrent 
contributions are included via the microscopic Riccati 
amplitudes a and b. 

For the derivation of the result ([TOjh selfconsistency has 
been assumed. However, even if the current-phase rela- 
tion used to evaluate Eq. (fTOj) has not been calculated 
selfconsistently, an upper bound for the change of the 
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FIG. 2: The volume V for the integration of the free energy 
encloses the Josephson junction (JJ) through which a total 
current / flows. Cross sections of the superconductor far from 
the junction are planes of constant gauge-invariant phase. 
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FIG. 3: (Color online) Three exemplary current-phase rela- 
tions (upper panel) and the corresponding variations of the 
free energy (lower panel). A: I = 7 c sin7, a normal or 0- 
JJ with 70 = 0; B: / = — 7 c sin7, a 7T-JJ with 70 = 7r; C: 
/ = — J c sin(27), a so-called cp-JJ with < 70 < 7r. 



free energy follows. Eq. (pTUj) is valid at arbitrary temper- 
ature as well as in the presence of an external magnetic 
field. 

To the best of our knowledge, the derivation of Eq. ([TO]) 
given in the present work is the first microscopic deriva- 
tion with general validity. Previous derivations were ei- 
ther based on thermodynamic reasoning and the appli- 
cation of the (second) Josephson relation dj/dt = 2eV/h 
with the voltage V across the junction^ or were re- 
stricted to tunneling junctions 26 . The derivation given 
in the present work does not depend on the actual re- 
alization of the JJ but is valid for arbitrary structures 
exhibiting a current-phase relation. 

By examination of Eq. (fTOj) , one finds that a zero of the 
l + ab current-phase relation /(70) = with ^(7) 1 7=70 

> 

corresponds to a stable local minimum of the free energy, 
and thus yields the intrinsic phase difference 70. If there 
exists no more than one nontrivial zero with < 70 < 7r, 
four cases can be distinguished: (1) 70 = corresponds to 
a normal JJ; (2) 70 = tt corresponds to a 7T-JJ. Finally, 
< 70 < 7r and I c > (I c < 0) corresponds to a so- 
called (^-J J 1Q i 27 with a positive (negative) critical current, 
where the critical current I c is defined as the absolute 
maximum of the current-phase relation. 

In Fig. [3l we sketch three exemplary current-phase re- 
lations and the corresponding variations of the free en- 
ergy. In the case of a normal JJ with 70 = 0, the curva- 
ture of the variation of the free energy at 7 = is positive, 
d 2 £(7)/d 7 2 | 7 = > (see curves A in Fig. [3]). For a ir JJ 
with 70 = 7r, however, the curvature of the variation of 
the free energy at 7 = is negative, d 2 £(j) /dj 2 \ 7 ^o < 
(see curves B in Fig. [3]). In the crossover regime with in- 
termediate intrinsic phase differences < 70 < 7r, higher 
harmonics dominate, but still the curvature of the varia- 
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tion of the free energy at 7 = is negative (see curves C in 
Fig. [3]). Accordingly, in the case of a normal or 0-JJ, the 
free energy of the junction firstly increases with increas- 
ing phase difference 7 > 0, whereas it firstly decreases in 
the case of tt and ip JJs. 



IV. SELFCONSISTENT SOLUTIONS 

In this section, we present full two-dimensional selfcon- 
sistent solutions for the microbridge geometry depicted 
in Fig. [TJ Therefore, we numerically calculate a selfcon- 
sistent solution of the selfconsistency equation for the 
pairing potential A(r, kp), Eq. (pQ). Based on this solu- 
tion for the pairing potential, we numerically solve the 
equation for the current density j(r), Eq. ([2]). 
For the selfconsistent calculations, we assume a cylindri- 
cal Fermi surface with the cylinder axis aligned perpen- 
dicular to the film plane. The geometry used for the cal- 
culations spreads over an area of about 12.5x12.5^0 with 
the coherence length £0 = Hvf / {^^oo{T = 0)). Specu- 
lar boundary conditions lead to n • j = with the sur- 
face normal n at all surfaces of the geometry. For the 
left and right end of the geometry depicted in Fig. [H 
periodic boundary conditions have been used. Selfcon- 
sistency automatically guarantees current conservation, 
V • j = 0. Details of the selfconsistent calculations have 
been published in a previous work on s-wave supercon- 
ducting microbridges 24 . 

In Fig. [H we present selfconsistent configurations for 
the amplitude and the phase of the pairing potential as 
well as for the corresponding current density. For this 
figure, an orientation angle of the d-wave of a = 7r/4 has 
been used because, in this case, the effect of the d-wave 
symmetry is most pronounced. The opening angle of the 
wedge has been chosen to be (3 = and the width of the 
microbridge w = 3.14 ^o- An intermediate temperature 
of T = 0.5 T c has been used. 

For the orientation angle a = 7r/4, the <i-wave symme- 
try leads to a suppression of the amplitude of the pairing 
potential at all surfaces of the rectangular geometry (see 
Fig. [4] (a)). The phase difference 7 across the microbridge 
which has been used for this figure corresponds to the 
critical current. At T = 0.5 T c , the phase of the pairing 
potential monotonically increases from —7/2 at the left 
boundary of the geometry to +7/2 at the right boundary 
(see Fig. 21(b)). The corresponding current distribution 
in Fig. [4] (c) exhibits contributions which flow along the 
gradient of the phase (from left to right, in positive direc- 
tion) as well as backflowing surface current a 18 i 19 . These 
backflowing surface currents are directly related to the d- 
wave symmetry and they are carried by Andreev bound 
states which exist at surfaces of d-wave superconductors 
with orientation angles a 7^ Q 15 i 16 i 17 . However, if one 
integrates the current density shown in Fig. [H (c) over 
a cross section of the geometry, a positive total current 
follows. 

In Fig. \5\ we present selfconsistent solutions corre- 
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FIG. 4: (Color online) Amplitude (a) and phase (b) of the 
pairing potential as well as the corresponding current density 
(c) at the critical current for a d-wave orientation angle of 
a = 7r/4, an opening angle of the wedge of /3 = and a width 
of w = 3.14 £0 at a temperature of T = 0.5 T c . The boundary 
geometry is indicated by the thick (red) lines. In (c), the 
shading is proportional to the current density. 



sponding to those in Fig. [H but for a lower temperature 
of T = 0.1 T c . From Fig. [5] (a), it is obvious that the am- 
plitude of the pairing potential hardly changes between 
T = 0.5 T c and T = 0.1 T c . However, the configuration of 
the phase of the pairing potential changes completely (see 
Fig. 03(b)). With decreasing temperature, the influence 
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FIG. 5: (Color online) Amplitude (a) and phase (b) of the 
pairing potential as well as the corresponding current density 
(c) at a temperature of T — 0.1 T c . Other parameters and 
presentation corresponding to Fig. \b\ 



of surface Andreev bound states strongly increases, which 
leads to dominating backflowing surface currents (Fig. [5] 
(c)). Accordingly, the phase exhibits a nonmonotonic 
variation with a phase shift (Fig. [5] (b)). At T = 0.1 T c , 
integrating the current density over a cross section of the 
geometry yields a negative total current. 

From the selfconsistent configurations presented in 
Figs. [H and it follows that the behavior of the mi- 
crobridge strongly depends on the temperature T. With 



decreasing temperature, the influence of surface Andreev 
bound states increases and, accordingly, the relative 
weight of the backflowing surface currents. An increas- 
ing importance of surface Andreev bound states with de- 
creasing temperature has also been found in a previous 
work on the influence of the surface Andreev bound states 
on the Bean-Livingston barrier 28 . The backflowing sur- 
face currents dominate at low temperatures T as well as 
for small widths w of the junction. This leads to negative 
values of the current-phase relations 29 and even to neg- 
ative critical currents 2 ^. As follows from Eq. ([T0|) which 
links the current-phase relation and the variation of the 
free energy, these negative currents are related to finite 
intrinsic phase differences < 70 < tt. The intrinsic 
phase difference 70 determines the state of the micro- 
bridge J J, which can be either a normal or 0-JJ, a cp-JJ 
or a 7T-JJ. 

The two-dimensional selfconsistent solutions shown in 
Figs. [H and [5] each are representative configurations for 
one set of the parameters T, a, (3 and w and for a fixed 
value of the phase difference 7. Fig. [4] for T = 0.5 T c 
corresponds to a normal or 0-JJ whereas Fig. [5] corre- 
sponds to a 7T-JJ. In order to obtain current-phase re- 
lations and, accordingly, intrinsic phase differences for 
diverse combinations of the relevant parameters, the full 
two-dimensional calculations have to be repeated many 
times. 

Finally, it should be noted that the selfconsistent con- 
figurations shown in Figs. [4] and [5] do not change sub- 
stantially if the orientation of the d-wave deviates from 
a = 7r/4. Small deviations do not lead to an abrupt 
disappearance of the backflowing surface currents. Sim- 
ilarly, opening angles of the wedge other than (3 = do 
not lead to an abrupt disappearance of the backflowing 
surface currents and do not substantially change the cur- 
rent configurations as presented in Figs. [Hand [H 



V. STEP MODEL 

In order to calculate complete phase diagrams of the 
microbridge, the selfconsistent calculations presented in 
the last section would have to be repeated for arbitrary 
combinations of the relevant parameters: temperature T, 
width of the junction w, orientation angle of the <i-wave 
a, opening angle of the wedge (3 and, finally, in order to 
obtain current-phase relations, phase differences across 
the junction in the range < 7 < tt. Unfortunately, 
the numerical costs of the full two-dimensional selfcon- 
sistent calculations inhibit the selfconsistent calculation 
of complete phase diagrams. Therefore, in the following, 
we employ a non-selfconsistent step model for the pairing 
potential A(r, kp) in order to calculate current-phase re- 
lations and the according variations of the free energy for 
the microbridge. Nevertheless, based on the microscopic 
derivation of the relation between the current-phase re- 
lation and the variation of the free energy, Eq. ([T0|h we 
know that a non-selfconsistent calculation provides an 
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FIG. 6: (Color online) Boundary geometry for the step model, 
Eq. (|12p . together with an exemplary quasiparticle trajectory 
r(s) and the Riccati amplitudes a(s), b(s) corresponding to 
Eqs. (H3l) . 



upper bound for the variation of the free energy. A de- 
tailed comparison of the full two-dimensional selfconsis- 
tent calculations and the step model will be discussed 
subsequently. 

We assume a cylindrical Fermi surface with the cylin- 
der axis aligned perpendicular to the film plane. Ac- 
cordingly, vp = ^(xcos# + ycos#) with x, y being 
unit vectors in the film plane and being the polar 
angle. Thus, in the case of d-wave pairing, Vk F ,k' = 
V cos(2# — 2a) cos(26 >/ — 2a) . The step model corresponds 
to an opening angle of the wedge of (3 = and assumes a 
step-like variation of the phase of the pairing potential, 
whereas its amplitude is taken to be constant: 

A L)jR (r,k F ) = A oc (T)cos(2^-2a)e^/ 2 (12) 

Here, the indices L,R label the left and right side of 
the junction and Aoo(T) is the temperature-dependent 
amplitude of the pairing potential in the bulk. The 
step model (|T2|) has to be solved taking into account 
the boundary geometry defining the microbridge JJ (see 
Fig. El). 

In order to find the current density ([TTj) at the cross sec- 
tion of the constriction, we solve the Riccati equations 
along trajectories r (s) = (xo = 0, yo) + s (cos sin 0) , see 
Fig. [U Introducing 



n(0) = V4 + |Aoo(T) cos(2#-2c0| 2 



■ ft(0) 



as well as rj(6) = 2Q(0)/(Hvf) and I = \yo/ sin0|, we find 
for < 6 < tt/2: 



a(s = 0) — cll{0) 
b(s = 0) = b R {9) 



(13) 



b R (27r-e)-b R (e) ^ 2Q(9) 



Results for tt/2 < 6 < 2tt follow accordingly. 
Based on Eqs. ([T3|) . the total current ([TT]) can be cal- 
culated by integrating the current density over the cross 
section of the microbridge J J, i.e. along the negative y- 
axis of the geometry depicted in Fig.[6l The integral over 
the cross section has to be taken from y = to the width 
of the junction at y = — w. 
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FIG. 7: (Color online) Phase diagram of the geometric 
Josephson junction shown in Fig. [T] The thick (red) line sep- 
arates regions of positive and negative critical current (above 
and below). For this figure, a = tt/4. 
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FIG. 8: (Color online) Intrinsic phase difference 70 corre- 
sponding to the phase diagram in Fig. [7] For this figure, 
a — 7r/4. 



VI. RESULTS AND DISCUSSION 

Based on the calculation of current-phase relations and 
critical currents I c from the step model ([T2|h we show in 
Fig.0the phase diagram of the geometric microbridge J J 
from Fig.[T]for a fixed d-wave orientation angle a = 7r/4. 
From the phase diagram, it follows that if the width of 
the junction is smaller than a critical value w c , the critical 
current is negative. For T — > 0, we find a value of about 
w c w 3.6^o- With increasing temperature, the critical 
width decreases to about w c « 0.78 ^0 near T c . Near T = 
T c , the current-phase relations assume the asymptotic 
forms / = J c sin7 for w > w c and / = — 7 c sin7 for 
w < w Cj respectively, and only the and the tt state 
occur. With decreasing temperature, higher harmonics of 
the current-phase relations become more important and 
the ip state appears in the vicinity of the 0-7T transition. 
At low temperatures, the f state extends to widths w 
much larger than the critical width w c which separates 
L < from L > 0. 
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FIG. 9: (Color online) Phase diagram for the variation of the 
d-wave orientation angle a. The thick (red) line separates 
regions of positive and negative critical current (above and 
below). For this figure, T = 0.5 T c . 

In the limit w — > 0, all trajectories traveling through 
the junction suffer a reflection. Accordingly, this situa- 
tion can be considered as a tt point contact, the comple- 
mentary configuration to a normal point contact^. The 
limit w — > implies I — > and Eqs. (|T3j) become partic- 
ularly simple: 

a(s = 0) = a L (0) b(s = 0) = b R (2ir - 0) (14) 

In this case, the current-phase relations of the <i-wave 
point contact are being reproduced, but with an intrin- 
sic phase shift of 70 = tt. From the phase diagram in 
Fig. [3 one finds that the tt point contact exists at all 
temperatures < T < T c . 

In Fig. [HI we show the intrinsic phase difference 70 
corresponding to the phase diagram in Fig. The dis- 
appearance of the (f state near T = T c becomes apparent 
as a discontinuous transition from 70 = to 70 = tt. 
At lower temperatures, however, a continuous transition 
arises. 

In Fig. [9j we plot the phase diagram for a fixed tem- 
perature of T = 0.5 T c , focusing on the variation of the 
orientation of the d-wave a. Starting from the ideal ori- 
entation a = 7r/4 for the occurrence of the tt state, we 
find that the critical width w c decreases when the d-wave 
is being rotated. However, small deviations from a = tt/4 
do not lead to an abrupt disappearance of the tt or the <p 
state which is important for the experimental realization. 

In Fig. [lOl we show current-phase relations and the 
corresponding variations of the free energy for <i-wave 
orientation angles a from a = to a = tt/ 4. The tt state 
is apparent for a close to 7r/4 since negative currents 
occur for all 7. Accordingly, close to a = 7r/4, the free 
energy decreases with increasing 7. With decreasing a, 
a transition to the state occurs, with positive currents 
and an increasing free energy for all 7. Close to the 
very transition, a ip region with an additional zero of the 
current-phase relation occurs. Since the gradient of the 
current-phase relation at this additional zero is positive, 




0.0 0.2 0.4 0.6 0.8 1.0 



y/n 

FIG. 10: (Color online) Current-phase relations (upper panel) 
and the corresponding variations of the free energy (lower 
panel) for angles a from a = to a = tv/4 (as indicated) in 
steps of 7r/100 for w = 0.1 ^o- The scale for the currents is 
given by I = 7veN(0)v F k B T c ^od. For this figure, T = 0.5 T c . 



it corresponds to a stable energetical minimum and an 
intermediate intrinsic phase difference < 70 < tt. 

In the present work, we use the step model ([T2|) in 
order to calculate phase diagrams of the geometric mi- 
crobridge Josephson junction depicted in Fig. [TJ From 
the microscopic derivation of the relation between the 
current-phase relation and the variation of the free en- 
ergy, we know that an upper bound for the free energy 
follows if a non-selfconsistent model for the pairing po- 
tential is employed. Based on a detailed comparison of 
current-phase relations from the step model and from 
full two-dimensional selfconsistent solutions 29 , we do not 
expect selfconsistency to qualitatively alter the general 
properties of the device as described here. We rather 
find that the critical width w c which marks the transi- 
tion to negative critical currents is underestimated in the 
scope of the non-selfconsistent step model. It should be 
noted that, according to the selfconsistent calculations, 
the opening angle of the wedge (3 hardly influences the 
current-phase relations. However, small opening angles 
of the wedge could possibly lead to an increased capaci- 
tive coupling of the electrodes. As long as the width w is 
in the range of several £0 or below, the step model proves 
to be a useful approximation. Since microscopic sur- 
face roughness does not suppress surface Andreev bound 
state o 18 i 31 , the reported intrinsic phase differences are 
expected to be robust features. 
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FIG. 11: (Color online) Dc superconducting quantum inter- 
ference device (dc SQUID) geometry consisting of two geo- 
metric microbridge Josephson junction according to Fig. [T] 
For a = 7r/4, the geometric microbridge Josephson junction 
labeled 1 can be in the 0, ip or ty state depending on its width 
and on the temperature whereas the Josephson junction la- 
beled 2 always is in the normal or stated. 



VII. CONCLUSION 

In the present work, we provide a microscopic deriva- 
tion of the relation between the variation of the free en- 
ergy and the current-phase relation with very general va- 
lidity. Based on a comparison with full two-dimensional 
selfconsistent solutions, we use a step model for the pair- 
ing potential in order to calculate current-phase relations 
and intrinsic phase differences for the geometric micro- 
bridge Josephson junction depicted in Fig. [TJ The cal- 
culation of intrinsic phase differences is used to access 



phase diagrams of the Josephson junction, justified by 
the relation to the variation of the free energy. 

From the phase diagrams, we conclude that the struc- 
ture sizes required for the experimental realization of 
the 7r and in particular of the cp state in cuprate high- 
temperature superconductor microbridges are within 
reach of modern fabrication technology. Because of the 
larger coherence length, electron-doped materials are es- 
pecially promising 32 . To test our predictions we suggest 
an interference experiment with a dc superconducting 
quantum interference device (dc SQUID) consisting of 
two microbridge JJs with ct\ = 7r/4 and a 2 = 0, re- 
spectively (see Fig. [TT]). Geometry- and temperature- 
dependent intrinsic phase differences according to the 
phase diagram will show up as shifts of the correspond- 
ing flux-dependent interference pattern. The experimen- 
tal verification of the intrinsic phase differences would at 
the same time imply a direct confirmation of the anoma- 
lous counterflowing quasiparticle surface currents which 
are a unique and intriguing fingerprint of d-wave pairing 
symmetry. 
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